arXiv: 1509.02780vl [hep-ph] 9 Sep 2015 


Prepared for submission to JHEP 


MS-TP-15-15 


Soft gluon resummation for associated ttH 
production at the LHC 


Anna Kulesza,“ Leszek Motyka,^ Tomasz Stebel^ and Vincent Theeuwes“ 

Institute for Theoretical Physics, WWU Munster, D-4-8149 Munster, Germany 
^Institute of Physics, Jagellonian University, S.Lojasiewicza 11, 30-348 Krakow, Poland 
E-mail: anna.kulesza@uni-muenster.de, leszekm@th.if.uj.edu.pl, 
tomasz.stebel@uj.edu.pl, vthee_01@uni-muenster.de 

Abstract: We perform resummation of soft gluon corrections to the total cross section 
for the process pp —>■ ttH. The resummation is carried out at next-to-leading-logarithmic 
(NLL) accuracy using the Mellin space technique, extending its application to the class of 
2—7-3 processes. We present an analytical result for the soft anomalous dimension for a 
hadronic production of two coloured massive particles in association with a colour singlet. 
We discuss the impact of resummation on the numerical prediction for the associated Higgs 
boson production with top quarks at the LHC. 
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1 Introduction 

Establishing the properties of the Higgs boson discovered at the LHC in 2012 [1, 2], in 
particular its couplings to the Standard Model (SM) particles, is one of the main tasks 
of the current LHC run. Since the SM Higgs boson couples to fermions proportionally to 
their masses, the top-Higgs Yukawa coupling is expected to be especially sensitive to the 
underlying physics. A direct way to probe the strength of the coupling without making any 
assumptions regarding its nature is provided by the measurement of Higgs production rates 
in the pp —)• tiH process. Although the production cross section is low and the collision 
energy and the luminosity available so far have not been sufficient enough to measure a Higgs 
signal in Run 1 [3-7], such a measurement in Run 2 is eagerly awaited. Correspondingly, 
precision predictions for the pp —)• tiH production process are of great importance and a 
lot of effort has been invested in the recent years to improve the theoretical accuracy. 

The next-to-leading-order (NLO) QCD, i.e. O^a^a) predictions are already known for 
some time [8, 9] and have been newly recalculated and matched to parton showers in [10- 
13]. As of late, the mixed QCD-weak corrections [14] and QCD-EW corrections [15] of 
are also available. Furthermore, the NLO QCD corrections to the hadronic tiH 
production with top and antitop quarks decaying into bottom quarks and leptons have been 
recently obtained [16]. Concurrently, new methods for a better measurement of the process 
have been proposed e.g. in [17] or in [18]. In general, for the LHC collision energies of Run 
2, the NLO QCD corrections are ~ 20%, whereas the size of the (electro)weak correction 
is more than ten times smaller. The scale uncertainty of the NLO QCD corrections is 
estimated to be ~ 10% [8, 9, 19]. While matching hxed-order predictions to parton showers 
pursued recently by many groups in such frameworks as aMC@NLO [10, 11, 20], POWHEG 
BOX [12, 13, 21] or SHERPA [22] allows for a more accurate description of hnal state 
characteristics, it does not change the predictions for the overall production rates. An 
improvement in the accuracy with which these rates are known can only be achieved by 
calculating higher order corrections. However, calculations of the next-to-next-to-leading- 
order corrections are currently technically out of reach. It is nevertheless interesting to 


- 1 - 




ask the question what is the size and the effect of certain classes of corrections of higher 
than NLO accuracy. In particular, we focus here on taking into account contributions from 
soft gluon emission to all orders in perturbation theory. The traditional (Mellin-space) 
resummation formalism which is applied in this type of calculations has been very well 
developed and copiously employed for description of the 2—7-2 type processes at the Born 
level. The universality of resummation concepts warrants their applications to scattering 
processes with many partons in the final state, as shown in a general analytical treatment 
developed for arbitrary number of partons [23, 24]. Recently, the soft gluon resummation 
technique in the soft collinear effective theory (SCET) framework was applied to pp —)• ttW^ 
[25]. So far, however, no calculations in the traditional resummation framework for processes 
involving 2 —3 scattering at the Born level have been performed. 

In this paper we take the first step in this direction by developing the Mellin-space 
threshold resummation formalism at the next-to-leading-logarithmic (NLL) accuracy for 
the case of 2 —)• 3 processes with two coloured massive particles in the final state. We 
then apply this formalism in order to estimate the impact of soft gluon corrections on the 
predictions for the total ttH production rate. In this particular case, the threshold region is 
reached when the square of the partonic center-of-mass (c.o.m.) energy, s, approaches M = 
2mt+mH, where mt is the top quark mass and mn is the Higgs boson mass. In the threshold 
region, the cross section receives enhancement in the form of logarithmic corrections in 
/3 = y^l — M^/s. The quantity (3 measures the distance from absolute production threshold 
and can be related to the maximal velocity of the tt system. Additionally, in the threshold 
region the virtual QCD corrections are also enhanced due to Coulomb-type interactions 
between the two final state top quarks which become large when the top quark velocity in 
the tt c.o.m. frame fiki —t 0 with fiki = y^l — Am^/ski and s^i = {pt + Pi)'^- However, the 
contributions to the total cross section from the threshold region are strongly suppressed 
by the /3^ factor originating from the massive three particle phase space. Nevertheless, 
one expects that the threshold corrections can still have a non-negligible impact on the 
predictions. 

The associated production of a Higgs boson with a tt pair involves four coloured partons 
at the Born level and as such is characterized by a non-trivial colour flow. The colour 
structure influences the contributions from wide-angle soft gluon emissions which have to 
be included at the NLL accuracy. The evolution of the colour exchange at NLL is governed 
by the one-loop soft anomalous dimension [23, 26-29]. Starting from four coloured partons 
in the process, the soft anomalous dimension is a matrix and is known for heavy-quark 
[26, 30], dijet [27-29] and supersymmetric particle production [31, 32], as well as for the 
general case of 2 —?■ n QCD processes [23, 24]. Here we adopt the calculations of the soft 
anomalous dimension for the case of 2 —>• 3 processes with two coloured massive particles 
in the final state. 
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2 Resummation for 2 —3 processes with two massive colored particles 
in the final state 


The resummation of soft gluon corrections to the total cross section Upp^ttH is performed 
in Mellin space, where the Mellin moments are taken w.r.t. the variable p = M'^/S. At the 
partonic level, the Mellin moments for the process ij —?• klB^ where i,j denote massless 
coloured partons, k, I two massive quarks and B a massive colour-singlet particle, is given 

by 


= dp p^ k-R) (2-1) 


with p = 1 — /3^. 

At LO, the ttH production receives contributions from the qq and gg channels. We 
analyze the colour structure of the underlying processes in the s-channel color bases, 
{c?} and {cf}, with c\ = c| = cf = 

Cg^ = • In this basis the soft anomalous dimension matrix becomes diagonal in 

the production threshold limit [26] and the the NLL resummed cross section in the A^-space 
has the form [26, 30] 


^(res) _y^;:^(0) 
/ 


^'ii^'klB.N—? ,^ii^klB.I.N Cij^klB,! ^V+1 


AJ A(“t) 
^N+l^ij^klB,I,N+l^ 


( 2 . 2 ) 


where we suppress explicit dependence on the scales. The index I in Eq. (2.2) distinguishes 
between contributions from different colour channels. The colour-channel-dependent contri¬ 
butions to the LO partonic cross sections in Mellin-moment space are denoted by / jv' 

The radiative factors A)y describe the effect of the soft gluon radiation collinear to the ini¬ 
tial state partons and are universal. Large-angle soft gluon emission is accounted for by 
the factors which depend on the partonic process under consideration and the 

colour conhguration of the participating particles. The expressions for the radiative factors 
in the MS factorisation scheme read (see e.g. [30]) 


1 _ 1 /■M2(1-z)2 j 2 

‘ -\Aiia,{q^)), 


lnA)v= / dz—^ -^ ‘ ^ '"2 

Jo 1 — z 


^^^ij^klB,I,N—J^ 1 — Z 


Dij^klB {1 — z)‘^)). 


The coefficients Aj, Dij^kiBj are power series in the coupling constant cts 


TT TT 


Dij^kW,! = ( — 
TT 




+ ... 


(2.3) 


The universal LL and NLL coefficients are well known [33, 34] and given by 

A^^ = Ci, Af = 1 Ci ((fl - Ca - |nj) with Cg = CA = 3, and C, = Cp = 4/3. 

The NLL coefficients Dij^kiB,! are obtained by taking the threshold limit s —)■ = 

{nik + rrii + thb)'^ oi the gauge-invariant soft anomalous dimension matrices In 

this limit = ^diag( 7 ^'^,...) and Dij^kiBj = 2Re(7y). The calculations of 
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apply the methods developed in the heavy quark heavy quark pair-production [26] to the 
process at hand, taking into account 2—7-3 kinematics, and yield 


pgq- 


MB 


TT 


CY{Lp,kl 

2^3 


1 ) 


— Q-j. 


5[(Ca ~ ‘^CF){Lfi^ki + 1) + C'aAs -|- (SCf — 3CA)ti3] 


, (2.4) 


with 


Ygg^klB 


Ois 

TT 


Tf? 0 O 3 

0 rii fns , 

_2ti3 %^f^3 rff _ 


(2.5) 


r?f = -CF(L^,ki + 1 ), 

rff = rff = ^((Ca - 2Cp}(L^,ki +1) + CaAs), 

where 


As = {Ti{mk) + T 2 {mi) + Ui{mi) U2{mk))l2, 

Os = {Ti{mk) + T2{mi) - Ui{mi) - U2{mk))/2, 


and 




Li3,kl 

Ti{m) 

Ui{m) 


+ l^ki fki- Pkl 


1 


2Kl3ki Ik + /3m * 


K = 


(Pi-Pkf, 


- (ln((m^ - tif/{rr?s)) -I + i-k) , 

^ (ln((m^ - Uif‘/{m^s)) - 1 + iTr) , 

\/l- {nik - rniY/sku Ski = {pk +Plf-, 

h = {Pj -Plf, Ui = {pi - pif, U2 = {pj - Pkf- 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


Eqs. (2.4), (2.5) reproduce the known results for heavy quark-antiquark (squark-antisquark) 
pair- production soft anomalous dimension [26, 31] in the limit ps —t 0. Also, our result 
for agrees with the result obtained in the SCET framework in [25]. It can be also 

explicitly seen that in the limit s —?• + the non-diagonal elements vanish and the 

diagonal elements give Dqq^kiBj = {0,-A^c}, Dgg^kiBj = {0,-A^c,-A^c}, which are the 
same coefficients as for the heavy-quark pair production Dij^kl- This confirms a simple 
physical intuition that the properties of the soft emission in the absolute threshold limit 
are only driven by the colour structure of the subprocesses and do not depend on the their 
kinematics. 

The coefficients 

_ 1 I I 

contain all non-logarithmic contributions to the NLO cross section taken in the threshold 
limit. More specifically, these consist of Coulomb corrections, A^-independent hard con¬ 
tributions from virtual corrections and A^-independent non-logarithmic contributions from 
soft emissions. Although formally the coefficients Cij^kiB,! begin to contribute at NNLL 
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accuracy, in our numerical studies of the pp —)• ttH process we consider both the case of 
Cij^kiB,! = 1) he. with the first-order corrections to the coefficients neglected, as well as the 
case with these corrections included. In the latter case we treat the Coulomb corrections 
and the hard contributions additively, i.e. 

^(1) _ ^(l,hard) ^(l,Coul) 

For k,l denoting massive quarks the Coulomb corrections are = C'Fvr^/(2/3fci) 

and = (Cp — C'A/2)7r^/(2/3fc;). The additive treatment is consistent with NLL 

resummation and matching to NLO. We note that in general Coulomb corrections can 
also be resummed. A combined resummation of Coulomb and soft corrections is, however, 
beyond the scope of this paper. 


3 Theoretical predictions for the pp —)■ ttH process at NLO+NLL accu¬ 
racy 


The resummation-improved NLO+NLL cross sections for thepp —?■ ttH process are obtained 
through matching the NLL resummed expressions with the full NLO cross sections 


a 


(NLO-l-NLL)^ 
h\h2^kl 


(p, pU M?,)ft. A) + ft. Pp) 


(res—exp) 


with 


a 


(res—exp) _ 

h\h2^klB ' 


~(res) / 2 2 \ - (res) / 2 2 \ i 

HjJkwAP' t^R) - HjJkiBMP^ Pf, Pr) \ 


\(NLO) 


(3.1) 


where is given in Eq. (2.2) and \(nlo) represents its perturbative expan¬ 

sion truncated at NLO. The moments of the parton distribution functions (pdf) fi/hix, Pp) 
are defined in the standard way f^p^ip'r) = Jq dx The inverse Mellin 

transform (3.1) is evaluated numerically using a contour C in the complex-Ai space accord¬ 
ing to the “Minimal Prescription” method developed in Ref. [35]. 

As mentioned in the previous section, the calculation of first-order contributions to the 
coefficients requires knowledge of the NLO real corrections in the threshold limit as 

well as virtual corrections. In our calculations we follow the methodology of [36, 37], where 
the case of two massive coloured particle in the final state was considered. We have explicitly 
checked that adding a massive colour singlet particle in the final state does not introduce any 
extra terms dependent on the mass of the added particle. Thus the Ai-space results for the 
pair-production process of two massive coloured particles are also applicable in our 2 —)• 3 
case. This way, the problem of calculating the j coefficients reduces to calculation of 

virtual corrections to the process. We extract them numerically using the publicly available 
POWHEG implementation of the ttH process [13], based on the calculations developed 
in [9]. The results were then cross-checked using the standalone MadLoop implementation in 
aMC@NLO [10]. Since the qg channel receives only colour-octet contributions, the extracted 
value contributing to is exact. In the gg channel, however, both the singlet 
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and octet production modes contribute. The implementation of the virtual corrections to 
gg —)• ttH in POWHEG and in aMC@NLO does not allow for their separate extraction 
in each colour channel. Instead, we extract the value which contributes to the coefficient 
^(i,hard) Qygj. gglour channels and use the same value to further calculate 

and In order to measure the size of the error introduced by this procedure, we 

then rescale this value by the ratios of the corresponding colour-channel dependent and 
colour averaged coefficients found for gg —)• tt in [38]. The scale dependence of the i 

can be fully deduced from renormalization group arguments, in the same way as for the full 
NLO result. We have checked that numerical results obtained with the procedure which 
we use to extract the values of the coefficients at go = gp = gp show the same scale 
dependence as expected from exact analytical expressions. 

In our phenomenological analysis we use nit = 173 GeV, mn = 125 GeV and choose 
the central scale gpg = gpfi = mt + mu/2, in accordance with [19]. The NLO cross section 
is calculated using the aMG@NLO code [20]. In the implementation of the resummation 
formula, Eq. 2.2, we numerically take a Mellin transform of the LO cross sections and the 
I coefficient terms which are both calculated in the x space. We perform the current 
analysis employing MMHT2014 [39] pdfs and use the corresponding values of ag. Beside 
presenting the full result including non-zero cj^j]^ttH I coefficients, we also show the results 
with = 1- 

We begin our numerical study by analysing the scale dependence of the resummed total 
cross section for pp —)• ttH at ^/S = 8 and 14 TeV, varying simultaneously the factorization 
and renormalization scales, gp and gp. As demonstrated in Fig. 1, adding the soft gluon 
corrections stabilizes the dependence on g = gp = gp of the NLO+NLL predictions with 
respect to NLO. As an example, the central values and the scale error at ^/S = 8 TeV 
changes from 132^g'®^ fb at NLO to Idl^^’^^ fb at NLO+NLL (with Ctj]^ttH i coefficients 
included) and correspondingly, from 613('(g^^" fb to fb at \/~S = 14 TeV. It is also 

clear from Fig. 1 that the coefficients strongly impact the predictions, especially at 

higher scales. 

In order to understand these effects better, in Fig. 2 we analyse the dependence on the 
factorization and renormalization scale separately for the case study of \/H = 14 TeV. We 
observe that the weak scale dependence present when the scales are varied simultaneously is 
a result of the cancellations between renormalization and factorization scale dependencies. 
A similar effect of the opposite behaviour of the total cross section under gp and gp 
variations was previously shown for the total cross section for the inclusive Higgs production 
in the gluon-fusion process [40]. The typical decrease of the cross section with increasing 
gp originates from running of ag. The behaviour under variation of the factorization scale, 
on the other hand, is related to the effect of scaling violation of pdfs at probed values of 
X. In this context, it is interesting to observe that the NLO+NLL predictions in Fig. 2 
show very little gp dependence around the central scale, in agreement with expectation of 
the factorization scale dependence in the resummed exponential and in the pdfs cancelling 
each other, here up to NLL. The relatively strong dependence on gp of the NLO+NLL 
predictions with non-zero C^j^ttH I then easily understood: the resummed expression 
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will take into account higher order scale dependent terms which involve both j and 

logarithms of N. These terms do not have their equivalent in the pdf evolution since the 
pdfs do not carry any process-specihc information. Correspondingly, they are not cancelled 
and can lead to strong effects if the coefficients j are numerically substantial. As 

these terms can only provide a part of the full scale dependence at higher orders, it is to 
be expected that their impact will be signihcantly modihed when NNLO corrections are 
known. 

Given the arguments above, we choose to estimate the theoretical uncertainty due to 
scale variation using the 7-point method, where the minimum and maximum values obtained 
with = (0.5, 0.5), (0.5,1), (1,0.5), (1,1), (1,2), (2,1), (2, 2) are considered. 

The effect of including NLL corrections is summarized in Table 1 for the LHC collision 
energy of 8, 13 and 14 TeV. The NLO+NLL predictions show a signihcant reduction of the 
scale uncertainty, compared to NLO results. The reduction of the positive and negative 
scale errors amounts to around 20-30% of the NLO error for = 13,14 TeV and to 
around 25-35% for y/S = 8 TeV. This general reduction trend is not sustained for the 
positive error after including the j coefficients. More specihcally, the negative error 

is further slightly reduced, while the positive error is increased. The origin of this increase 
can be traced back to the substantial dependence on /Xi? of the resummed predictions with 
non-zero coefficients, manifesting itself at larger scales. However, even after the 

redistribution of the error between the positive and negative parts, the overall size of the 
scale error, corresponding to the size of the error bar, is reduced after resummation by 
around 7% at 8 TeV and 10 (13)% at 13 (14) TeV with respect to the NLO uncertainties. 
The scale error of the predictions is still a few times larger than the pdf error, cf. Table 1. 
For simplicity, the pdf error shown in Table 1 is calculated for the NLO predictions, however 
adding the soft gluon correction can only minimally influence the value of the pdf error. 

As expected on the basis of large phase-space suppression in the threshold regime, the 
predictions for total cross section at NLO+NLL are only moderately increased by 2-3% 
w.r.t. the full NLO result. Introducing the coefficients j leads to an increase in 

the iL-factor of up to 6-7%, indicating the importance of constant terms in the threshold 
limit. Since the impact of soft corrections is bigger for processes taking place closer to 
threshold the iL-factor gets slightly higher for smaller collider energies. We also check 
the impact of our approximated treatment of keeping parts of and 

coefficients coming from the virtual corrections equal to the colour channel averaged value, 
by rescaling at fip = the averaged coefficient with ratios 

taken from [38]. The procedure is motivated by obvious similarities between the colour 
structures of the pp —)• tt and pp —)• ttH cross sections considered at threshold. We hud 
that such rescaling of the hadronic ttH cross section leads to a 3 per mille effect at 14 
TeV, or a 5% effect on the correction itself. Therefore we do not expect that the exact 
knowledge of the i and C^gg^^m g coefficients will have a signihcant impact on the 

hadronic NLO+NLL predictions. However, we stress that because of the large phase-space 
suppression in the threshold regime the resummed results, while systematically taking into 
account a well dehned class of correction, should not be used to estimate the size of the 
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Vs [TeV] 

NLO [fb] 

NLO+NLL 

NLO+NLL with C 

pdf error 



Value [fb] 

K-factor 

Value [fb] 

K-factor 


8 

1 00+3.9% 

1 qc:+3.0% 

1.03 

141+7.7% 

1.07 

+3.0% 

-2.7% 

13 

p-nf;+5.9% 
0 UD 94 % 


1.02 

c q 7 + 8 . 2 % 
^^^-5.5% 

1.06 

+2.3% 

-2.3% 

14 

q+6.2% 

01^-94% 

e9c:+4.6% 

1.02 

DDU _5 7 % 

1.06 

+2.3% 

- 2 . 2 % 


Table 1. NLO+NLL and NLO total cross sections for pp — >■ ttH for various LHC collision energies. 
The error ranges given together with the NLO and NLO+NLL results indicate the scale uncertainty. 


NNLO total cross section, by e.g. methods of expansion of the resummed exponential. 
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Figure 1. Scale dependence of the LO, NLO and NLO+NLL cross sections at '/S = 8 and 
\/S = 14 TeV LHC collision energy. The results are obtained while simultaneously varying pp and 

PR, p = pp= PR. 


4 Summary 

We have investigated the impact of the soft gluon emission effects on the total cross section 
for the process pp —)• ttH at the LHC. The resummation of soft gluon emission has been 
performed using the Mellin-moment resummation technique at the NLO+NLL accuracy. 
To the best of our knowledge, this is the first application of this method to a 2 —s- 3 process. 
Supplementing the NLO predictions with NLL corrections results in moderate modifications 
of the overall size of the total rates. The size of these modifications, as well as the size of 
the theoretical error due to scale variation is strongly influenced by the inclusion of the 
first-order hard matching coefficients into the resummation framework. The overall size 
of the theoretical scale error becomes smaller after resummation, albeit the reduction is 
relatively modest when the non-zero first-order hard matching coefficients are considered. 
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_1_1_1_ 
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Figure 2. Factorization and renormalization scale dependence of the LO, NLO and NLO+NLL 
cross sections at \/~S = 14 TeV LHC collision energy. The results are obtained with = fj,Q for fj,p 
variation and and fip = fiQ for fip variation. 
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